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A one-dimensional electron system at quarter-filling has been examined by applying the renor- 
malization group method to a bosonized model with on-site ([/) and nearest-neighbor {V) repul- 
sive interactions. By evaluating both normal scattering and Umklapp scattering perturbatively, 
we obtain a phase diagram in which a metallic state with a 2fcF spin density wave (fcr is the 
Fermi wave number) moves into an insulating state with charge disproportionation of a 4fcp 
charge density wave with an increase in both U and V . The effect of the next-nearest-neighbor 
repulsion is also discussed. 
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A one-dimensional (1-D) electron system at quarter- 
filling is a basic model for understanding th£ electronic 
properties of quasi 1-D organic conductors The 1-D 
electron system with only on-site Coulomb repulsive in- 
teraction, U (1-D Hubbard model), does not exhibit a 
metal-insulator (M-I) transition as a function of J7, and 
is metallic at any filling, except half-filling.B The most 
dominant state is given by 2k-p spin density wave {2k-p- 
SDW), where k-p is the Fermi wave number. The pres- 
ence of the long-range repulsive interaction is expected to 
enrich the phase at quarter-filling. In fact, numerical di- 
agonalization for the modal jvith both U and the nearest- 
neighbor interaction, y,oD'B'B' shows that the insulating 
phase appears for a large strength of both U and V and 
that the superconductivity becomes the most dominant 
fiuctuation for a large V and small U . On the other 
hand, several proDepties have been elucidated within the 
mean-field theory.Q'&B With increasing V ^ a transition 
occurs from a pure 2fcF-SDW state to a coexistent staJje 
of 2fcF-SDW and Ak^ charge density wave (4A:f-CDW)E) 
Such a coexistence is maintained for the charge ordering 
in (DLJ)CNQI)2Ag observed by the ^^C-NMR measure- 
ment .US' The transition has been examined by evaluating 
the commensurabilitv energy corresponding to the 8fcF- 
Umklapp scattering.^ The next-nearest-neighbor repul- 
sion results in a coexistence of 2fcF-SDW and purely elec- 
tronic 2/cf-CDWJj^ and it has been proposed to be the 
origin of the coexistencej-observed by the X-ray experi- 
ment on (TMTSF)2PF6.0 

At quarter-filling, the 8fcF-Umklapp scattering is cru- 
cial to obtain the commensurabiUty energy. Although 
the above mean-field results may show common features, 
we need to calculate the Umklapp scattering in the 1-D 
system by taking into account quantum fluctuation. The 
existence of commensiixahilit\!_energy of such high order 
has been pointed out£3'E3olij-' However, to the best of 
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our knowledge, there are no studies for the phase dia- 
gram on the plane of U and V , which is calculated by 
using the analytical expression of the commensurability 
energy. 

In the present paper, a 1-D system with repulsive 
interactions at quarter-filling is investigated using the 
bosonization method and the rcnormalization groiip 
(RG) theory. Based on the commensurability energyjES 
a phase diagram is derived on the plane of U and V . 
The relevance of the present results to the observation 
in (DI-DCNQI)2Ag sah,EHl' and the effect of the next- 
nearest-neighbor interaction are discussed. 

We consider a 1-D extended Hubbard model given by 
the Haniiltonian, Ti ~ Tio + TYint, 



Ho = —t 



JO- 



h.c.j - ^y^ri-j^g 



(1) 



Ka 



Hint — 



^J,(7^^-f 1,(7' 



a a' Ki~A4 ^ ^ 
X 5Ki+K2-K3-Ki,Ga^K^,a°''^K2,<j''^K3,a'aKi,c! , (2) 

where t and /i denote the transfer energy and chemi- 
cal potential, respectively, and eK — —2t cos Ka with 
lattice constant a and —n/a < K < tt/u. The quan- 



operator of the electron at the j-th site with spin cr, 
rij^a — o^ja^h<^^ ^ ~ 0, ±27r/a, and Nl is the num- 
ber of the lattice. For the later convenience, we divide 
the one-particle states as dk,-,a = aK.a for —■n/a < 
K < -7r/(2a), Ck-,a = aK,a for -7r/(2a) < K < 0, 
Cfe,+.<T = o-K,a for < K < 7r/(2a), and dk,+,a — o-K.a for 
7r/(2a) < K < n/a, k being the deviation of the wave 



tityaj g.(= 1/\/Nl X^k^ '^""'at-,) denotes the creation 
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number from ±kp for Ck,±,a, and iSfep for dk.±.a-iS^ 
In terms of Ck,p,cr and dk.p.a, 'Hq is written as Ho — 

and Hint is rewritten as Tiint — X]i=o^int,i, where i de- 
notes the nmiiber of dk^p.a and/or d\ ^ the respective 
interactions. 

When the one-particle states only near the Fermi wave 
number, zbfcp, are taken into account, the Hamiltonian 
is reduced to 



pka 



pkk qaa 



|(5i±^<t-<t' + gi\\Saa')cl.^g.p^^cl.,_^_p^^,Ck',p,a'Ck,-p,a- 
+ {92±^cr-a' + g2\\'^aa')cl^qp^^cl,_^_p^^,Ck',-p,<T'Ck,p,a 
+ igi±Scr-a' + 9i\\^<Ta')cl+q^p^„cl,_^^p^^,Ck'^p,a'Ck,p,cr^ , 

(3) 

where L — N^a and the energy dispersion is linearized at 
ibfep as e±fcp+fc — ji ~ ivpk with vp = ^/2ta. Coupling 
constants are given by g';l^ = = Ua/2 + Va, g^y = 
Va, gjj^ = Ua/2, and g°^^ = 0. We note that 
eq. (3) does not show the M-I transition since there is no 
Umklapp scattering. I— . 

According to a theoretical suggestion by SchulzH and 
calculation of the mean-field theory,® Umklapp scatter- 
ing in the case of quarter-filling appears through the 
one-particle states with ztS/cp, which are connected to 
those of the Fermi point by ±2fcF. We take into ac- 
count the states near iSfep systematically by the follow- 
ing procedure in order to derive the effective Hamilto- 
nian, which consists of the one-particle states only near 
i/cp. Representing Ck^p^a and dk^p^a- in terms of Grass- 
mann algebra, the partition function of eqs.(l) and (2), 
Z, is given by Z = / 2?Mfe,p,^dfe,p,^]2?[cfc^p^^Cfc,p,^]e-^, 

where S = So[cl^p^^Ck,p,a] + Soldi j,^^dk^p^a]+J2t=o •^'mt,*, 
and Siat,i denotes the action corresponding to 7iint,i- 
After integrating dk,p,a, the partition function, Z, 
is written as Z = J T^lcl ^ ^Ck,p,cr]e~^'ff , where 
Z'l = /2?[d^,p,.4,p,.]e-^°[<->'=— 1, and S^ff = 

So[clp,aCk.,p,a] - ln(exp(-^1'^(3S'i„t,i))^ ((••■)<i is the 
average by So[dl p^^dk^p,cr])- By calculating Seff pertur- 
batively, the effective Hamiltonian is derived where the 
diagram for the vertex includes only the green functions 
of dk,p,a- From the perturbation up to the second order 
of U and/or V , the correction to the normal processes. 



ill(-L) 



in eq.(3) is calculated as 



Ag2\\ = -2Di [Vaf , 



Ag,, = -2D,(^) -2dJ^ 



Va 



A.g4|| 
Ag4± 



-2D2{VaY , 
2D,\^-Va 



(4) 
(5) 
(6) 
(7) 
(8) 



where L>i = (87rto)-i /J^^^ d2/(siny + 1/^2)-^ = 
(87^^a)"l^/21n(^/2 + 1) ~ 1.25/(87rta) and D2 = 
>(47ria)-i /J'^^dj/(sin?/+cos?/+\/2)-i = {%'!ita)-^2V2/ {^2+ 
1) ^ 1.17/(87rta). The 8fcp-Umklapp scattering, which 
conies from the third order expansion, is obtained as 

Wi/4 = ^^(f/a-4Fa) 



{Uaf 
8i2 



^ j Ax{ipl,^^_p^„f'ipl_^ip_p,. 



where the first term coixesponds to the conventional 
commensurability energyllj) andV'p.a = ^/y^J2k^^''^'^k,p,(T- 
In the third order, there are other contributions leading 
to the correction to the normal processes, which are dis- 
regarded in the following due to eqs.(4)-(8) being large 
enough for the present choice of parameters. _ 

Here we utilize the bosonization methocO by in- 
troducing phase variables for the charge (spin) fluc- 



tuation, 9p and 



and (j)cr), which are defined 



by V) = W(9i)e-"l^l/^-'^^E;op.('^)4+,.p..Cfe,p,. 
and(/.p(^^) = W(9i)e;"l'l/'-'';Efep.Mp4+,,p,<.Cfc,p. 
respectively. The quantity a ^ is the ultraviolet cut- 
off. Those fields satisfy, [9^{x), (f>^' {x')] — i7rsgn(a; — 
x')S^^i {v — p or a). The electron operator is ex- 
pressed as ?/'p,o- — (27ra)^^/2exp{i(p/2)(6'p-|-p0p-|-cr0cr + 
(jp4>a)} exp(i7rE:p^) with Sp+ = p/2 Y.p' Np>+ and Sp_ = 

vl^Yjp' ^v'- + Ep' ^P'+ where Np„ = Jdx^pl^i^p^. In 
terms of these phase variables, the Hamiltonian is ex- 
pressed as H = Hp 



Tin 



-P. 

47r 



Ha + H' , where 



gl/4 

2(7ra)2 



dx cos 40 



Ha — 



H' = 



9i± 

2{'KOif 



1 



dx{^{d^daf 



Ka{d^CPaf 



dx cos 26a , 
dx cos 4:9 „ cos 29n 



(10) 

(11) 
(12) 



Here = By j Ay and Vy — VY\JByAy {v — p or 
cr), where Ap^a) = 1 + {.94|| + (-).94i + g2\\ + (-).92± - 
9i\\}/{t^vf) and Sp(^) = 1 + {g^ + {~)g4± - 52|| - 
{+)92± + 9i\\}/{t^vf) with = 5° + Ag^i^^y In 

eqs.(10)-(12), the coefficients of the nonlinear terms are 
respectively given by gi± = y^j^ + Agi± = (?7a/2){l — 
ADi{Ua/2 - Va)}, gi/4 = {Ua)^{Ua - 4ya)/(27ri>p)2 x 
(a/a)2 and 5^/4 = (Uaf /{2TTVFf/2x {a/ af. For F = 0, 
the quantity gi/^-is essentially thCpSame as that re- 
ported by Schulal^) and Yonemitsu.O A contribution 
from the foiLrth order of U has also been suggested by 
GiamarchiJiil^ The order parameters of the 2fcp-SDW, 



correlation Hjttccts in a une-uimensionai i^uarter-t luea Hjiectron system witn nepuisive interactions 



2fcF-CDW and 4fcF-CDW are given by 
02fep-SDW = 5]ae-'2?"=F-^t_^^_^^^ 

pa 

oc s\n{2kYX + Op) sin da , (13) 
02fep-CDW = ^e-'2p'=F-^t^^^_^^^ 

pa 

oc cos{2k^x 6 p) COS 6a- , (14) 

04fep-CDW = 5]e-'4f'=F-^t^^vi-V'-p.-V'-p.+ 
p 

oc cos(4fcFa; + 29p) . (15) 

We investigate the possible states in the hmit of low 
energy. Since the scaling dimension of 7i' given by 
2 — 8Kp — 2Kcr is smaller than that of the other non- 
linear terms, the term may be safely neglected for de- 
termining the phase diagram. In this case, the Hamilto- 
nian is divided into the charge part and spin part. For 
the spin part, as long as gi± > 0, the quantities, gi± 
and i^CT tend toward and 1, respectively, for the low 
energy limit, due to SU(2) symmetry, and the excita- 
tion is gapless. Then the long range correlation func- 
tions are given as (sin 6*0.(2;) sin0o.(O))pj ln^/^(a;) and 
(cos6'o(a;)cos6'<,(0)) - x"Mn"^/^(x)B The low energy 
property for the charge part is determined by the follow- 
ing RG equations, 

j^Kpil) = -^Gl^mlil) , (16) 

^Gi/4(0 = [2-8ifp(/)]Gi/4(0 , (17) 

with the initial conditions given by Kp{Q) = Kp and 
^1/4(0) = gi/4^ll2'n:Vp). Here I = ln(Q;'/Q!) with the new 
length scale a' larger than a. The solution of the above 
RG equations is obtained by 

Gl/S)-^l{'^Kp{l))-2\nKp{l)^ const , (18) 

with a decrease of Kp with increasing I. 

The RG flows are_shown in Fig.l, where the value of 
a ~ 2a /tt is usedO In the regions (I) and (II), the 
nonlinear term becomes irrelevant, indicating a metal- 
lic state. The correlation functions with long distance 
are given by {sa\Op{x)sa\Op{Q)) ~ (cos 0p(a;) cos 0p(O)) ~ 
^-K,{oo) g^^^ (cos20p(a;)cos26lp(O)) ~ x-^^'-^^\ There- 
fore the dominant state exhibits a crosspter between 2fcF- 
SDW and 4fcF-CDW at Kp{oo) = 1/30) where Kp{oo) 
is determined by Gl^^{Qi) - l/{2Kp{<S)) - 2\iYKp{Q) = 
— l/(2ifp(oo)) — 21ni4rp(oo). The phase boundary deter- 
mined by Kp{oo) — 1/3 is shown by the dashed curve. 
Both regions, (III) and (IV) are insulating states where 
Gi/4 — > —00 and 6p is locked as 0,7r/2 mod tt in (III), 
and Gi/4 00 and 9p — ±7r/4 mod tt in (IV). The or- 
dered state of 4A:f-CDW leading to charge ordering is 
realized in (III) and 6p — {Op — tt/2) corresponds to 
the state in which the charge is rich at the even (odd) 
sites. The closed circle, open circle and closed square 
correspond to the initial condition of y/t = 0, 2 and 
6, respectively, with U /t = 5, showing that the metallic 
state with 2fcF-SDW in (I) moves to the insulating state 




Fig. 1. Scaling flows ( dotted curves ) derived from eqs. (16) and 
(17) on the plane of Kp and G1/4, where the solid and dashed 
curves denote their boundaries. The parameters in regions (I) 
and (II) lead to weak coupling, where Kp{oo) = 1/3 on the 
dashed curve. A strong coupling regime is found in both regions, 
(III) and (IV). The closed circle, open circle and closed square 
correspond to the initial values of V/t = 0, 2 and 6, respectively, 
with the flxed U/t = 5. 
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Fig. 2. The phase diagram on the plane of V/t and U/t, where 
regions (I), (II) and (III) separated by boundaries (solid and 
dashed curves) correspond to those in Fig. 1. The metallic 
states with 2fcF-SDW and 4fcF-CDW are obtained in regions (I) 
and (II), respectively, while region (III) is the insulating state 
with Ak-p charge ordering. The dotted curves in the metallic 
state denote Kp{oo). 



with the 4fcF-CDW ordering with increasing V. 

The phase diagram on the plane of U/t and V/t is 
shown in Fig. 2. The insulating phase appears for a large 
strength of U and V. In addition, for large V and small 
U, Kp{oo) approaches unity indicating that the super- 
conducting fluctuation is enhanced in this region. Figure 
2 is qualitatively the sama as^that derived from the nu- 
merical diagonalization.Q'lj'ElEP In the insulating phase, 
which corresponds to the region (III) in Fig.l, the 4fcF- 
CDW, i.e., charge disproportionation, is realized. The 
state is consistent with the mean-field resultfl) Thus, the 
insulating state obtained by the numerical diagonaliza- 
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tion is due to the formation of charge disproportionation, 
which originates in both U and V interactions. Here, we 
discuss the effects of the nonhnear term, cos 49 p cos 29a^ 
of eq.(12), which may give minor changes on the phase 
boundary. In the insulating region, the term seems to 
give the correction of gi± as gi± — > gi±+g[/4^ (cosiOp) /2 
with {cos 46 p) > 0. Such a couphng constant gives rise to 
the spin gap. However, the procedure breaks the SU (2) 
symmetry obviously. Therefore, it is expected that the 
third order correction of K^, which has been discarded 
in the present treatment, would restore the symmetry, 
and the spin excitation remains gapless in the insulating 
region. 

Within the present analysis, the insulating state ex- 
hibits the fluctuation of 2A:f-SDW, while 4A:f-CDW ex- 
ists as the true long range order. When three dimension- 
ality is introduced by the interchain coupling, the fluc- 
tuation of 2fcF-SDW becomes the true order. However, 
the characteristic temperature of 2/cf-SDW is expected 
to be smaller than that of 4A:f-CDW. This conclusion 
may be relevant to the following experimental observa- 
tion in (DI-DCNQI)2Ag. The magnetic order of 2fcF- 
SDW is observed at 5.5 K, whereas the resistivity shciffis 
insulating behavior at least below room temper ature.LZP 
In addition, charge disproportionation is observed below 
220 K.y) 

Finally, the effects of the next-nearest-neighbor repul- 
sion, V2, are briefly discussed. Up to the second or- 
der perturbation in the present formulation, the prop- 
erties of the charge fluctuation do not show a qualitative 
change. On the other hand, the coefficient of the non- 
linear term of the spin fluctuation, gi±, is changed as 
{Ua/2-V2a){l-iDi{Ua/2-Va + V2a)}. Since the non- 
linear term with positive (negative) gi± favors 2fcF-SDW 
(2fcF-CDW), it is obvious that the repulsive next-nearest- 
neighbor interaction stabilizes 2fcF-CDW. In fact, the 
spin excitation becomes gapped and 2fcF-CDW is real- 
ized when U/2 < V2, independent of the strength of V 
as seen in the mean- field result Thus the quantity V2 
stabihzes (suppresses) 2fcF-CDW (2fcF-SDW). However, 
the coexistence of two kinds of 2/cf density waves, which 
is predicted by the mean-field theoryja) cannot be un- 
derstood from such a calculation, which needs further 
exploration. 

In conclusion, we investigated both normal scattering 
and Umklapp scattering analytically for one-dimensional 
quarter-filled electron systems with on-site and nearest- 
neighbor interactions. We have shown that the increase 
of both U and V leads to the insulating state with charge 
disproportionation and that the next-nearest-neighbor 
interaction suppresses (enhances) the fluctuation of 2fcF- 
SDW (2fcF-CDW). 
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